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Abstract. Theory is presented for the distributions of local process intensity and 
local average pore dimensions in random fibrous materials. For complete partitioning 
of the network into contiguous square zones, the variance of local process intensity 
is shown to be proportional to the mean process intensity and inversely proportional 
to the zone size. The coefficient of variation of local average pore area is shown 
to be approximately double that of the local average pore diameter with both 
properties being inversely proportional to the square root of zone size and mean process 
intensity. The results have relevance to heterogenous near-planar fibrous materials 
including paper, nonwoven textiles, nanofibrous composites and electrospun polymer 
fibre networks. 
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1. Introduction 

The global average pore size of thin, i.e. near-planar, heterogeneous fibrous materials 
and its distribution have been widely studied using statistical geometry and simulation. 
Typically the context of these studies has been materials with widespread application 
in society and industry such as paper [1,2], nonwoven textiles [3], and fibrous filter 
media [4,5]. Interest in the structural characteristics of fibrous materials has increased 
in recent years as researchers seek to develop materials for future applications such 
as carbon nanotube 'buckypaper' [6,7], nanofibrous composites [8,9] and electrospun 
fibrous networks for application as scaffolds in tissue engineering [10, 11]. 

In a seminal paper, Miles [12] provides several properties of the polygons generated 
by the stochastic division of a plane by a Poisson process of straight lines with infinite 
length. A graphical representation of such a process is given in Figure la. The process 
intensity is characterized by the expected line length per unit area, f . The following 
results of Miles are utilized in the present study: 

• The expected area of polygons is 

a = - 2 - (1) 

• The distribution of diameters of the largest circle inscribed within 
polygons is exponential with mean 

d=\ (2) 

• a and d are independent of the width of lines, for any probability density 
of line widths. 

Miles showed also that the expected number of sides per polygon is 4 and that the 
fraction of triangles is -P(3) = (2 — ^-) ph 0.355; Tanner [13] obtained the fraction of 
quadrilaterals, -P(4) ps 0.381. The fractions of polygons with more than 4 sides are not 
known analytically, but have been obtained by Monte Carlo methods [14, 15]. We note 
that very similar results were obtained by Piekaar and Clarenburg [4] from simulations 
of networks of fibres with finite length. Crain and Miles [14] observed the probability of 
n-sided polygons was well approximated by the Poisson variable (n — 3) with mean 
(n — 3) = 1. This result was used by Dodson and Sampson [16] to approximate 
the probability density functions for polygon areas and perimeters, assuming regular 
polygons. The coefficient of variation of polygon area was close to 2 and that of polygon 
perimeter was y/3/2; despite the assumption of regularity, these are close to the analytic 
results of Miles [12]. 

Now, each intersection between lines represents one of the vertices of four polygons 
and for a Poisson process of infinite lines, Miles [12] gives the expected number of 
intersections per unit area as f 2 /n, i.e. I /a. For the case of networks lines with finite 
length and constant width (Figure lb), the expected number of intersections per unit 
area was derived by Kallmes and Corte [17] for infinite networks; in practical contexts, 
this can be considered satisfied of samples are large relative to the dimensions of voids. 




Figure 1. Graphical representations of Poisson line processes in a unit square with 
intensity f = 200. a) infinite lines; b) straight lines with length 0.1; c) curved lines 
with length and radius of curvature 0.1. 



Recently it has been shown [18] that Kallmes and Corte's expression for finite lines is 
the same as that of Miles for infinite lines. Further, Berhan et al. [19] showed the same 
dependency by considering the probability of intersection of segments of finite length 
curved fibres {e.g. Figure lc); validity was confirmed via simulation studies for networks 
of fibres with sinusoidal curvature with differing sinusoidal frequency. Accordingly, we 
may be confident that theory describing a planar Poisson process of infinite length lines 
provides an appropriate reference to probe the void structure of networks of finite and 
curved fibres, greatly simplifying the analysis. 

Statistical geometric models for the polygon area and pore diameter distribution 
suggest that these are described by gamma distributions [20-22] and this has been 
confirmed through simulation [5] and by measurements on paper [23,24], nonwoven 
textiles [25] and electrospun polymer nanofibre networks [26]. Whereas these models 
and measurements characterize the global average pore dimensions and its distribution, 
preliminary inspection of Figure 1 reveals significant local differences in pore dimensions 
corresponding to the underlying variability in process intensity. Quantification of local 
variability in pore dimensions is the focus of this study. 

The variance of local process intensity, and hence that of local mass density 
was derived for a planar Poisson process of rectangles by Dodson [27] and has 
since been widely employed as the standard reference in quantifying variability in 
paper [28]. Relationships between network uniformity and fracture behaviour have 
been demonstrated for this material [29,30]. Similarly, flow and filtration efficiency in 
fibrous filters exhibit variability between regions [31,32]. In a theoretical treatment, 
Chatterjee [33] considered the local average elastic modulus of nanofibrous composites 
and demonstrated a strong sensitivity to the pore size of the constituent fibre network. 
Holzmeister et al. varied local pore size in electrospun polymer scaffolds for tissue 
engineering by generating networks of fibres with bimodal distributions of diameter 
and remarked that the tendency of cells to proliferate on these scaffolds was influenced 
by local void dimensions [34]. 

Here we consider the distribution of local process intensity for a Poisson process 
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of infinite lines. From this we obtain approximate probability density functions for the 
distribution of local averages of pore area and pore diameter in planar stochastic fibrous 
materials. 



2. Variance of local process intensity 

We seek the variance of local process intensity for a planar Poisson line process of mean 
intensity, f , partitioned into contiguous square zones of side, x. We denote the local 
process intensity within such zones, r. 

Before proceeding, it is helpful to obtain estimates of the likely range of f 
encountered in real materials. If the expected mass per unit area, or 'areal density', of a 
fibre network is, j3 (kg m -2 ), and the linear density of the constituent fibres is 5 (kg m" 1 ), 
then the expected process intensity, f = (3/5. For materials such as paper and nonwoven 
textiles, with mass per unit area between 20 and 100 gm~ 2 we expect f to be of 
order 10 2 mm" 1 ; for electrospun networks of polymer fibres with mass per unit area 
around 10 gm" 2 and fibre diameter around 1 /im we expect f to be of order 10 3 mm _1 . 

From Coleman [35] the probability density of the length, li, of random secants in a 
unit square is 

^ if < h < 1 

m = \ ^m~^ lfl <^v/2 (3 ) 

otherwise. 

with mean U = 4 (1 - V2 + 3 log(l + V2)) /(3 tt) « 0.946. 

The expected total length of lines in a square zone of side x is 

L = fx 2 . (4) 

We denote the length of secants in a square of side x, l x , with expected length, l x = li x. 
It follows that the expected number of secants in a square of side x is 

i x ti 

If the number of lines contained in a given square of side x is N, then the total line 
length in that square is 

N 

such that the local process intensity is 

-J ■ m 

It follows that the variance of local process intensity is given by 

= ^ • (8) 
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where the subscript x is included to denote that the variance depends on the zone size. 
Note also that cf 2 {L) = x 2 <7 2 (Li), where L x is the total line length in a unit square 
containing N secants. Thus, 

= ^ ■ (9) 



X 



Now, Li is a random variable obtained as the sum of the lengths of N independent 



and identically distributed secants in a unit square, i.e. L\ = ^2f =0 h,i, where is a 
continuous random variable with probability density given by Equation (3) and N takes 
integer values. The mean and variance of Li are given by [36] 

Li =Nh (10) 
a 2 (L 1 )=Na 2 (h)+l 1 2 a 2 (N) . (11) 

We assume that TV" is a Poisson random variable, such that cr 2 (N) = N and we have 

a 2 (Z 1 )=iv(a 2 (/ 1 )+I 1 2 ) (12) 

Now, a 2 ^) = P 1 -1 1 2 . So 

a 2 {Li) = NT 2 . (13) 
We obtain l\ as the second moment of the probability density given by Equation (3): 



V2 



3 



l 2 = llmdh^- . (14) 
Jo 71 

For completeness, we note that o- 2 (l\) = (3/n) — l\ fa 0.0593. 

Substituting Equations (13), (14) and (5) in Equation (9) yields our final expression 

for the variance of local process intensity for square zones of side x: 



3 f 



nli x 
f 

x 



(16) 



2.1. Distribution of local process intensity 



From the Central Limit Theorem, we expect the distribution of local process intensity, r 
to be well approximated by a Gaussian distribution if the expected number of secants 
in a square zone of side x is sufficiently large. For low intensity processes and at small x 
we anticipate that the distribution of local process intensity will exhibit a positive skew 
as a consequence of the underlying Poisson process for N. 

Derivation of the distribution of r has proved intractable, so here we estimate the 
skewness of the distribution by considering a Poisson process of secants in a unit square. 
The influence of changing f or x is therefore captured entirely by varying the expected 
number of secants in the unit square, N . From Equations (5) and (15), the mean and 



Spatial variability of fibrous materials 



6 



0.4 



•3 0.3 



£ 0.2 



2 "d 



5IU 



0.1 



I (Ki 



200, 



500^0, 

Ann, 



p 



0.1 0.2 0.3 

Skewness, yi(x), from Equation (If 



0.4 



Figure 2. Skewness of simulation data plotted against that calculated using 
Equation (18). Data labels give mean number of lines in a unit square; broken line has 
unit gradient. 



variance are r = Nl% and o~ 2 (t) = 3N/tt. Noting that l\ ~ l\ ~ 1 and given that 
< h < y/2, we assume that If w 1 such that the third central moment, /i3, of the 
probability density of r is approximated by that of the underlying Poisson process, i. e. 
/i 3 (r) m ^s(N) = N. The skewness of the distribution of r is therefore approximated 
by 

7i(t0 = (17) 

3 

7T2 



3V3ViV 

To test the validity of Equation (18) and to probe the distribution of r further, a Monte- 
Carlo simulation of random lines in a unit square was carried out using Mathematica [37]. 
The number of lines in a unit square was generated as a Poisson random variable with 
mean N. For each random line, the code generated a pair of points drawn from a 
standard uniform distribution to represent a coordinate through which the line passed. 
This coordinate was associated with an angle drawn from a uniform distribution on 
the interval [— 7r/2, 7r/2). The equation of the line was generated from the coordinate 
and the angle; its points of intersection with the perimeter of the unit square were 
calculated, allowing the length of the line to be computed. For each N, 10 6 unit squares 
containing N lines were simulated and the local process intensity computed as the total 
length of the lines in each square. As a check, the distribution of 10 6 random lines 
was compared with the probability density given by Equation (3); the mean agreed 
within 0.003 % and the variance agreed within 0.11 %. 
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Figure 3. Histograms showing the frequency of local process intensity, r, arising from 
Monte-Carlo simulations of random lines in a unit square. Heights of bars have been 
scaled to give unit area under the curve. Solid line represents PDF of fitted skew- 
normal distribution; broken line represents PDF of Gaussian distribution with same 
mean and variance as the data. 



The skewness of the local process intensity, 71 (t), of the simulation data arising 
from different input N is plotted against that obtained using Equation (18) in Figure 2. 
Data labels represent the value of N input to the simulation and the broken line has 
unit gradient. A linear regression on the data has gradient 0.9994 with coefficient of 
determination, r 2 = 0.9997. 

Histograms of the local process intensity arising from the simulations are plotted in 
Figure 3 for N = 10, 20 and 50. To approximate the probability density of the data, the 
heights of the bars are given by the frequency divided by the bin width. The solid lines 
represent the probability densities of skew-normal distributions fitted by a least-squares 
method to the cumulative data. The probability density of the skew-normal distribution 
is [38] 



erfc 



a (m— t) 



9{r) = * ^ ' , (19) 

'2 TT S 



where erfc(£) is the complementary error function. The mean, variance and skewness 
are given by 

\r2as 

r = m + 



TV y/1 + a 2 



ir (1 + or 



71(70 



3 



\/2(4-7r)a 
2^ 



7T + 7T 



respectively. 

As anticipated from the skewness values plotted in Figure 2, the distributions 
are increasingly well approximated by a Gaussian probability density as N increases. 
When is greater than about 50, the skewness is negligible and the Gaussian can be 
assumed to describe the distribution of local process intensity well. 



Spatial variability of fibrous materials 



8 




Local average pore area, a (fim 2 ) Local average pore area, a (\irn-) 



Figure 4. Probability density function for local average pore area a as given by 
Equation (22). Left: influence of process intensity f at scale of inspection x = 1 mm; 
right: influence of scale of inspection x when process intensity f = 200 mm -1 . 



2.2. Distribution of local average pore area 

We proceed assuming the distribution of local process intensity to be Gaussian with 
mean f and variance f/x (cf. Equation (16)), i.e. 



/ X x(t-t) 2 

9 ^ = V^ eT ' (20) 

We note that for fx > 20, g{r) dr < 4 x 10" 6 ; accordingly, truncation of the 
distribution such that < r < oo is unnecessary for practical purposes. 
From Equation (1), we expect the local average pore area to be 

2=4 • (21) 

Inevitably, when a network is partitioned into contiguous square zones, some polygons 
intersect the perimeter of the zone. The expected number of polygons intersecting the 
perimeter of a square zone of side x can be approximated as n perim = Ax/d = 4xf 
and the expected number of polygons in the square is n area = x 2 /a = x 2 f 2 /n. 
So, the expected fraction of polygons intersecting the perimeter of the square is 
approximately n per i m /n ai . ea = 47r/(f x). When f = 100 mm" 1 and x = 1 mm this 
fraction is 0.125 and n area > 3000 so it is reasonable to assume that Equation (21) 
provides a good measure of the local average polygon area. 

The probability density of local average pore area is obtained by a simple variable 
transform of Equation (20): 



p(a) 



da 



x (22) 



\a 3 f 



The probability density given by Equation (22) is plotted for a range of process 
intensities, f and scales of inspection, x in Figure 4. As anticipated, the distribution 
exhibits a positive skew and narrows with increasing f and x. A consequence of 
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neglecting skewness and using the approximation given by Equation (18) for the variance 
of local average process intensity is that the expected local average pore area obtained as 
a = J °° ap{a) da is slightly greater than n/f 2 though for f > 200 mm -1 and x > 1 mm 
the error is less than 2 %. The variance of local average pore area, is given by 



a 2 (a) 



{a — a) 2 p(a) da 



(23) 



It has not been possible to obtain a closed form solution to this integral though an 
analytic estimate can be obtained through consideration of Equations (5) and (21). 
From Equation (5) we expect the local average process intensity to be 

_ Nh 



x 

Substituting in Equation (21), we obtain 
irx 2 1 

a = 



l{ N 2 



such that 



a 2 x {a) 



71 X 



J 2 



a 2 x {l/N 2 ) . 



(24) 



(25) 



(26) 



Now, N is a Poisson variable and since P(N = 0) > 0, a 2 (l/N 2 ) is undefined. 
Typically, we expect N to be sufficiently large that P(N = 0) is negligible such that 
< 1/N < 1. A convenient approximation to the discrete Poisson probability function 
is the probability density of a Gamma distributed continuous random variable with 
variance equal to the mean. This probability density is given by 

-N 



T(N) 

such that the probability density of v — 1/ iV 2 is given by 

diV 



(27) 



r{v) 



dv 



(l+N/2) 



2T(N) 

The distribution has mean 9 — 1/ ((N — 1)(N — 2)) and variance 

AN- 10 



(28) 



<j\v) = a 2 (l/N 2 ) 



(N 
4 

iV5 



l) 2 (A^-2) 2 (A^-3)(A^-4) 
for TV > 20 



(29) 



Such that 
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Figure 5. Coefficient of variation of local average pore area CV x (a) plotted against 
mean process intensity f . 



Substituting for N from Equation (5) yields 




(31) 



Since a = 7r/V 2 , it follows that the coefficient of variation of local average pore area is 
approximated by 



We observe that the influence of zone size and process intensity is coupled such that the 
coefficient of variation depends only on the dimensionless product, xf. The coefficient 
of variation of local average pore area, as calculated via numerical integration of 
Equation (23) is plotted against mean process intensity in Figure 5. The solid lines 
represent the approximation given by Equation (32). We note that Schweers and 
Loffler [31] report a coefficient of variation of local flow velocity through a porous 
nonwoven filter of 0.3 at the 0.5 mm scale, consistent with the values of CV{a) plotted 
in Figure 5. 

2.3. Distribution of local average pore diameter 

Although the voids generated by random fiber processes are irregular convex polygons, 
it is often convenient to characterize their dimensions by an equivalent diameter, rather 
than by area. A good candidate for such a measure is the equivalent diameter determined 
from the hydraulic radius and defined as the ratio of the area of a polygon to its 
perimeter. Despite the established utility of this measure, to calculate it for our system 
we would require knowledge of the joint probability density of polygon perimeter and 
area, which is unknown. Two alternative measures have been employed previously: the 
diameter of the largest circle that can be inscribed within a polygon [5, 12] and the 
diameter of a circle with the same area as a polygon [1,20]. The expected diameter 




2 



(32) 
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Figure 6. Probability density function for local average pore diameter d as given by 
Equation (34). Left: influence of process intensity f at scale of inspection x = 1 mm; 
right: influence of scale of inspection x when process intensity f = 200 mm -1 . 



of inscribed circles is given by Equation (2), and from Equation (1) it follows that 
the expected diameter of circles with the same area as a polygons is d eq = 2/f. So, 
in both cases we see that the expected pore diameter is inversely proportional to the 
process intensity f , and we expect this dependency to hold for the local averages such 
that d = d eq /2 = 1/r. Accordingly, we proceed to derive the probability density for the 
local averages of diameters of inscribed circles and note that this can be readily scaled 
to give the probability density of equivalent diameters; for convenience, we will term 
the diameter of an inscribed circle the 'pore diameter'. 

The probability density of local average pore diameter is obtained by a variable 
transform of Equation (20): 



q(d) 



dr 



dd 



1 x 



g{l/d) (33) 

c (1-df) 2 

2«f»r . (34) 



d 2 V 2 7T f 

Following our earlier treatment, we estimate the variance of local average pore 
diameter. 



1 x 



(35) 



^ 2 M = X j I al{l/N) (36) 
n 

Again approximating the Poisson distribution for N by a gamma distribution with 
variance equal to the mean, we obtain 

> ( 37 ) 
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Figure 7. Coefficient of variation of local average pore diameter CV x (d) plotted 
against mean process intensity f . 



such that 



x 2 



„l(d) * (38) 



X T° 



and 



CV x (d)^\l— (40) 

XT 

' (41) 



Note that our estimate of the coefficient of variation of local average pore diameter 
given by Equation (41) applies for both pore inscribed pore diameters and the diameters 
of circles with the same area as pores. Comparing Equations (32) and (41) we observe 
that our approximation for the coefficient of variation of local average pore area is double 
that of the local average pore diameter, so again depends only on the dimensionless 
product fx. The coefficient of variation of local average pore diameter as calculated 
via numerical integration to obtain the moments of Equation (34) is plotted against 
mean process intensity in Figure 7. The solid lines represent the approximation given 
by Equation (41). 

In the introduction, we noted that the coefficients of variation of polygon perimeter 
and diameter were v3/2 and approximately 2, respectively, and thus are independent of 
process intensity. To a first approximation we expect the 'diameter' of convex polygons 
to be proportional to their perimeter, such that the coefficient of variation of diameters is 
constant also J. Here we have seen that the coefficients of variation of the local averages 
of pore diameter and pore area are not constant, and both exhibit the same strong 
dependence on the mean process intensity and the scale of inspection. 

| There is good experimental evidence that this is the case [1,39-41] 
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3. Conclusions 

Established theory for the global average polygon dimensions arising from planar Poisson 
line processes has been extended to give the distributions of local average polygon 
dimensions. Assuming that the number of lines passing through contiguous square 
regions of the network has a Poisson distribution, the variance of local average process 
intensity has been shown to be proportional to the mean process intensity and inversely 
proportional to the side of the square regions. Through an approximation, the skewness 
of the distribution has been shown to be inversely proportional to the square root of the 
expected number of lines passing through each region; this result has been confirmed 
by Monte Carlo simulation. 

For typical values of process intensity observed in heterogeneous fibrous materials, 
the skewness of the distribution of local process intensity rapidly approaches zero and 
the distribution can be assumed to be Gaussian. On this basis, the distribution of 
local average polygon area and that of local average pore diameter have been derived, 
along with approximate expressions for their coefficients of variation. We find that the 
coefficient of variation of pore diameter is inversely proportional to the square root of 
the product of process intensity and zone size; the coefficient of variation of pore area 
is double that of pore diameter. 
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